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Abstract 



J> ' The first law of black hole mechanics (in the form derived by Wald), is expressed in terms 

in ^ of integrals over surfaces, at the horizon and spatial infinity, of a stationary, axisymmetric black 

' hole, in a diffeomorphism invariant Lagrangian theory of gravity. The original statement of the 

, first law given by Bardeen, Carter and Hawking for an Einstein-perfect fluid system contained, in 

addition, volume integrals of the fluid fields, over a spacelike slice stretching between these two 
\C , surfaces. One would expect that Wald's methods, applied to a Lagrangian Einstein-perfect fluid 

' formulation, would convert these terms to surface integrals. However, because the fields appearing 

, in the Lagrangian of a gravitating perfect fluid are typically nonstationary, (even in a stationary 

<sJ black hole-perfect fluid spacetime) a direct application of these methods generally yields restricted 

5^ , results. We therefore first approach the problem of incorporating general nonstationary matter 

fields into Wald's analysis, and derive a first law-like relation for an arbitrary Lagrangian metric 
theory of gravity coupled to arbitrary Lagrangian matter fields, requiring only that the metric 
■ field be stationary. This relation includes a volume integral of matter fields over a spacelike slice 

rS , between the black hole horizon and spatial infinity, and reduces to the first law originally derived 

' by Bardeen, Carter and Hawking when the theory is general relativity coupled to a perfect fluid. 

We then turn to consider a specific Lagrangian formulation for an isentropic perfect fluid given 
by Carter, and directly apply Wald's analysis, assuming that both the metric and fluid fields are 
stationary and axisymmetric in the black hole spacetime. The first law we derive contains only 
surface integrals at the black hole horizon and spatial infinity, but the assumptions of stationarity 
and axisymmetry of the fluid fields make this relation much more restrictive in its allowed fiuid con- 
figurations and perturbations than that given by Bardeen, Carter and Hawking. In the Appendix, 
we use the symplectic structure of the Einstein-perfect fluid system to derive a conserved current 
for perturbations of this system: this current reduces to one derived ab initio for this system by 
Chandrasekhar and Ferrari. 



1 Introduction 



The first law of black hole mechanics as stated by Bardeen, Carter and Hawking Q relates small 
changes in the mass of a stationary, axisymmetric black hole to small changes in its horizon surface 
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area, angular momentum and the properties of a stationary perfect fluid that might surround it: one 
first fixes a stationary axisymmetric Einstein-perfect fluid black hole solution with stationary killing 
field (with asymptotically unit norm) and axial killing field (with closed orbits). One then 
defines b to be an infinitessimal perturbation to a nearby stationary axisymmetric solution; then the 
first law in |]^ is 

m = ^bA^^HbJH- [ IJ-'\v\5Nabc+ I ^6Jabc+ [ T\v\dSabc, (1) 

where the spacetime is characterised by an ADM mass, M, and the black hole by its horizon surface 
area. A, surface gravity, k, angular velocity, ^Ih, and angular momentum Jh (measured at the horizon). 
The fields associated to the perfect fluid are its four velocity, (which here is taken to be of the form 
jja _ ^Qy'l^l^ where v"" = (^'^ + ^lip*^, for some (generally non-constant) J7), the chemical potential iJ,', 
the temperature T, stress-energy T"^, and number and entropy densities n and S. The three-forms 
Nabc = nll'^eabcd, Jabc = T'^e^p'^edabc, and Sabc = SV^edabc represent the fluid number density, angular 
momentum density and entropy density on a spacelike 3-surface, S, that has boundaries at the black 
hole horizon and the two-sphere at spatial infinity. We have also set eabcd to be the canonical volume 
element on spacetime. 

Considerable effort has been spent on weakening the assumptions made in (|l|) on the background fields 
and their perturbations. For instance, consider an arbitrary diffeomorphism invariant Lagrangian 
theory with both metric and matter fields, and let the theory possess stationary, axisymmetric black 
hole solutions, which are asymptotically flat, and have a bifurcate killing horizon (for an explanation 
of these terms see [Q, ^). Then it was shown ||2|, ^, providing the metric and matter fields appearing 
in the Lagrangian were stationary and axisymmetric in the black hole background, that there existed 
a first law of black hole mechanics in a form only involving surface integrals on the sphere at spatial 
infinity and the bifurcation sphere of the black hole horizon. Namely, given the Lagrangian for the 
theory, one could algorithmically define integrals £ and J over the sphere at spatial infinity, and S 
over the bifurcation sphere, satisfying the following identity: 

5S = ^5S + VtH5J. (2) 
zvr 

(Here 5 denotes a perturbation from the background black hole solution to any nearby solution.) The 
quantity £ was interpreted as the canonical energy of the black hole system, J as the canonical angular 
momentum and S as the black hole entropy. 

We might therefore expect that the volume integrals in (||) involving the fluid can be converted to 
surface integrals in the form (^), by choosing a suitable variational form for the Einstein-perfect fluid 
system and using the methods of @. In fact, we are unable to reproduce the first law (|^) in a form 
only containing surface integrals, using these methods; the difficulty is that at least one of the fields 
appearing in each of the Lagrangian formulations for a perfect fluid (that we are aware of) is generally 
non-stationary, even when the fluid four velocity, number density, entropy, and functions of these fields 
(which we refer to collectively as the physical fields), are stationary. Since the methods of |^ require 
that all fields appearing in the Lagrangian (which we refer to henceforth as the dynamical fields) are 
stationary and axisymmetric in the black hole background, the allowed background solutions for the 
perfect fluid in the resulting first law are restricted. 

This paper gives two results in response to this problem: we first relax all explicit symmetry assump- 
tions on matter fields appearing in the Lagrangian, and find the consequence for the first law given in 
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1^] . We also attempt to generate a first law of the form (^) by a careful choice of an existing Lagrangian 
formulation for gravity coupled to a perfect fluid, directly using the methods of 

In section (2) we consider an arbitrary Lagrangian theory of gravity coupled to arbitrary matter fields, 
assuming only that the metric is stationary and axisymmetric in the black hole background, but 
making no such assumptions about the matter dynamical fields. We then modify the methods of 
to generate a perturbative relation, but instead of attempting to express the matter contribution to 
the first law (|2|) via surface integrals, we leave it instead as a volume integral over a hypersurface, S, 
joining the bifurcation sphere to the sphere at spatial infinity. In restricted cases (which we explain 
later) we can motivate an independent measurement of the "vacuum" black hole mass, Mg. In these 
cases we can also define quantities which resemble the "vacuum" black hole entropy, Sg, and angular 
momentum, Jgjf, and having done so our perturbative relation takes the form 

SMg = ^6Sg + nnSJgj,, + • e Saab -5{e - T- 0, (3) 

where T"^ is the stress-energy of the matter fields. We will see that this relation defines a black hole 
entropy, Sg, which is in general not the black hole entropy defined in [Q]: however, in special cases the 
interpretation of Sg as black hole entropy can be appropriate (for instance, as we show in section (4), 
this relation reduces to (|l|) when the gravitational theory is chosen to be general relativity, and the 
matter source is chosen to be a perfect fluid). Our result differs from a similar relation presented by 
Schutz and Sorkin [^, in that they conjectured, but did not explicitly include the black hole entropy 
and angular momentum boundary terms, and so did not explicitly generalise the full form of (|^). In 
addition, as we shall explain, the definition of our "Noether current" (involved in the intermediate 
calculations) is both less ambiguous than that presented by Schutz and Sorkin and more general 
than the definition given by Sorkin The range of theories in which our methods are well defined 
is therefore larger than those addressed by their methods. 

In section (3) we define a gravitating perfect fluid and review some variational principles for it: Schutz's 
"velocity-potential" formulation Q, which uses the dynamical fields ((/>, a, (3, 9, a) to define the product 
of the (physical) specific inertial mass and four velocity - iJ,Ua = + aVaP + 6'VaO', and Carter's 
more recent "axionic vorticity" formulation |jl^ for an isentropic perfect fluid, which uses a dynamical 
field bab to define the number current Nate (given in (||)) via Nate = 'i^[aKc]^ ^i^d the dynamical fields 
to define the fluid vorticity via 2V[a/uC4] = 2V[aX^V6]X • 

In section (4) we present two forms of the first law for the Einstein-perfect fluid system. The first 
form is derived from the relation (^) and is the same as (||), with the exception that 6 is now allowed 
to be a perturbation from the (stationary axisymmetric) background to an arbitrary nearby solution. 
(Note that this form of the first law contains volume integrals.) It is of also interest to know if we 
can construct any form of the first law with perfect fluids only involving surface integrals; in fact, by 
directly applying the methods of Q for a metric theory of gravity coupled to a perfect fluid described 
using Carter's variational principle, (with the potential bab for Nabc, and for '^ab) we can derive a 
first law of the form 

6M + fl^6 / bgr-fJ-o^S / bgr = -^SA - nnSJH + / - / Xgr, (4) 

where M is the ADM mass, A is the black hole surface area, Jh is the black hole angular momentum 
appearing in (|l|), Xgr is the two-form 2^Pbp[g[S{fiUrj) — Vr]X~Sx'^ + ^r]X'^Sx~]j and we have written 
S°° and TC for the sphere at spatial infinity and the bifurcation sphere, respectively. We will see that 
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this first law is more restrictive than (|l|), but it is the only non-trivial rule of the type involving a 
perfect fluid that we can currently construct. 



In the Appendix we evaluate the symplectic form of the Einstein-perfect fluid system, using the 
variational formulation given by Schutz |||] for the perfect fluid. The symplectic form is dual to a 
generally conserved current, quadratic in the field perturbations |11|. We find (in parallel with Burnett 
and Wald's calculation for the Einstein-Maxwell system [^]) that this conserved current reduces to a 
current previously derived ab initio by Chandrasekhar and Ferrari [^] for the polar perturbations of 
a static axisymmetric black hole. 



2 A perturbative relation for black hole mechanics with non-stationary 
matter fields 



In this section we give a perturbative relation that resembles the first law of black hole mechanics, 
for an arbitrary theory of gravity with a diffeomorphism invariant Lagrangian. We assume the the- 
ory possesses black hole solutions in which the metric is stationary and axisymmetric, but place no 
restrictions on the other fields appearing in the Lagrangian (we refer to these fields collectively as the 
dynamical fields). The motivation for this is, as we have indicated, that variational formulations for 
gravitating Einstein-perfect fluid systems have fluid dynamical fields which are nonstationary even 
when the fluid's physical fields (the four-velocity, number density and entropy) are stationary and 
axisymmetric. We first make some necessary definitions related to the the symplectic structure of a 
diffeomorphism invariant Lagrangian theory. These are explained in detail in ; here we merely state 
(and, in one case, refine) the relevant definitions and results. In the following we often use bold face 
type to denote differential forms on spacetime, suppressing their indices when convenient. 



2.1 Some Preliminaries 



All theories we consider arise from a Lagrangian, which is taken to be a diffeomorphism invariant 
four-form on spacetime, dependent on the metric, gab, and some arbitrary set of matter fields, ip. (We 
collectively refer to all the dynamical fields by <j).) By this we mean that the Lagrangian has the 
functional dependence 

L = L{gab, Rated, VRabcd, {^Y Rabcd, 1p,ViP, . . . , (V) V), (5) 

(here multiple derivatives appearing in the above expression are assumed to be symmetrised - see 
for further discussion about this dependence). In particular we require that every field appearing in 
the Lagrangian give rise to an equation of motion (there are no "background" fields). The variation 
of the Lagrangian defines these equations, E = 0, along with the symplectic potential 0, by 

dL = B6(l) + d@{(j),S(l)). (6) 

(Here @{(p,6(j)) is a linear differential operator in the field variations Sep. Because the Lagrangian is 
only defined up to the addition of an exact form, L — > L -|- dfj,, the symplectic potential is only defined 
up to the following terms: (0, 5(/>) @{(j),6(j)) + dY{(p,6(j)) + 5fj,{(f)), where Y and fj, are covariant 
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forms with the same type of functional dependence as and L, respectively. These ambiguities were 
discussed in Q.) 

Now fix a smooth vector field, on spacetime. Then the Noether current J[^] associated to is a 
three-form defined by 

J[e]^0(0,£^</')-^L (7) 

where the centred dot denotes contraction of the vector into the first index of the form. This Noether 
current can be seen to obey the following identity: 

dJ[e] = -BC^cj), (8) 

which we now use to further elucidate its structure. (Although they appear in a different context, the 
calculations below have the same fiavour as those in the Appendix of [^.) 

Lemma 1: Fix L to be the Lagrangian of a diffeomorphism invariant theory of gravity and matter 
fields, with equation of motion E = as given in (^). Without loss of generality, label each dynamical 
field (p hy i, and give each field Ui upper, and di lower indices: also label the equations of motion for 
each field similarly, so that the equation of motion term in (^) becomes 

E50 = 6£;^,,^...,^;--^>5,^,^--^"»,^...,^^. (9) 

Then for any smooth field ^" there exists a two- form, Q[^], called the Noether charge associated to 
(which is local in the dynamical fields and ^'^), such that the Noether current J[^], defined in (^), can 
be written 

J[^] = -(e ■ E ■ ■ + dQl^], (10) 

where we define the three-form 

i 



Proof: 

For clarity, we first consider the case where the metric is the only dynamical field: (j) gab- Then 
setting the metric field equations E"'' = eE"'^, Eq.(|8|) reads 

dJ[e] = -2eEfVa^b 

= -2eVa{Efib) + 2eVa{Ef)ib- (12) 

Therefore setting {e- Eg- Cjabc = edabcEg'^Ce, we have 

d(J[e] +2e-Eg-0 = 2eVa{Ef%, (13) 

which shows that the right side of (^) is both linear in and exact for all The results of Q now 
imply that the right side must vanish identically - and so V aE""^ = 0. This in turn implies that the 
left side of ( |l3|) must be an identically closed three-form, which (using the results of [^] again) implies 
the existence of a two- form, Q[^], local in the dynamical fields and such that 

m+'^^-Eg-i = dCt[i]. (14) 
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We define Q[^], the Noether charge associated to as any two form which is local in the dynamical 
fields and and satisfies this relation. 

We can also perform this analysis for L with the general dependence (^). With the labels for each 
field and its equation of motion given in @, the first equation in (|l^) becomes 

d3[(] = -eE^0..,....„;'^-'^'''A'/'.'^-'-a,...a,,, (15) 



which, through a similar manipulation to (12) leads to the structure for J[^] and the definition of the 
Noether charge, Q[C], in ([l0|). □ 

The Noether charge was defined in Q only when E = 0, via J[^] = dQ[.^]. This left open the definition 
of Q[^] when E 7^ 0. In the Appendix of |5|, however, it was shown that Q[^] could be defined when 
E 7^ 0, such that there existed forms Cq with J[^] — dQ[i^] = Ca^,"", and where the Ca vanished when 
E = 0. At that time it was not known whether Q[^] was uniquely defined this way, nor was the 
explicit form of Ca specified. We have given this explicit form in (11). Moreover, the above analysis 



uniquely defines the Noether charge via (|1C|), without imposing the field equations, up to the following 
ambiguities (which were discussed in detail in Q): The ambiguity in described after Eq. (P) means 
that J[^] is only defined up to the following terms: J[^] — > J[^] + d(Y((/), £^0) — ^ • fi), and so the 
ambiguity in Q[^] is Q[^] — > Q[^] + Y{cj},C^(j)) — ^ • /x. These ambiguities will not affect the results 
stated in the following sections. 



We now define the symplectic current, u((j),Si<j),624>), (a three-form on spacetime) by: 

u;((A, SicP, ^20) = S2@{<P, SicP) - 6i&{<P, 62(t>). (16) 
Note that a; is a function of an unperturbed set of fields, (p, and is bilinear and skew in pairs of 



variations {di(p,524>)- It can be shown (see [11|) that this three-form is closed when is a solution 
of the field equations and 5i(p and 524> are solutions of the linearised equations of motion (In the 
Appendix we examine this closed form - it is dual to a conserved vector field, which we evaluate for 
perturbations of an Einstein-perfect fluid system). Moreover, if we let be a smooth vector field, set 
6i(j) = C^(f) and let (520 = 50 be a variation to a nearby solution (with 55f^ = 0), then uj{(j),C^<j),6(j)) 
can be shown [Q] to be exact: 

c^(0, £^0, (50) = d[5Q[^] - e • 0(0, <50)]. (17) 



Now fix a black hole spacetime with a stationary and axisymmetric metric, for the theory given by 
the Lagrangian in (|5|); let the stationary killing field with unit norm at spatial infinity be and the 
axial killing field (with closed orbits) be (p"". Let the black hole have a bifurcate killing horizon, with 
bifurcation sphere 7i, and let it be asymptotically flat, with the two-sphere at spatial infinity S°°. 
Let S be a three-surface with these two boundaries, and set (50 to be an arbitrary perturbation of 
the background which satisfies the linearised equations. Then the first law of black hole mechanics as 
stated in Q is an interpretation of the identity 

/ ^(0,^50,(50)= / 5Q[^]-e-0(0,<50)- / (5Q[e]-^0(0,50) (18) 

(which arises from integrating (|l7|) over S). When Lie derives all the dynamical fields in the 
background, the left side of ( p!8|) vanishes, and one is left with a relation between surface integrals 
on the boundaries of S, which can be shown to be of the form (^). In section (4) we present an 
explicit Lagrangian for the Einstein-perfect fluid system, and, assuming that all dynamical fields are 
stationary and axisymmetric, compute the surface terms arising from this Lagrangian. 
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2.2 The perturbative identity 



Having stated these necessary definitions we turn to construct our perturbative identity. We start by 
decomposing the Lagrangian L into a part Lp, depending on the metric, gab, (which is assumed to be 
stationary and axisymmetric in the black hole background), and a part L^, dependent on both the 
metric and a set of matter fields, ip, (on which we place no restrictions): 

L = 'Lg{gab, Rabcd, ^Rabcd, ■ ■ ■ , (^YRabcd) 

+L™(V', VV', . . . , {Vfij, gab, Rabcd, V Rabcd, {^Y Rabcd)- (19) 

Since this breakup only requires that be independent of any matter fields, it is very non-unique, 
and in general we have no method of controlling the ambiguity 

+ 

Lm — > — A. (20) 

where A = \{gab, Rabcd, Rabcd, • • • , {"^y Rabcd)- 

The variation of the Lagrangian yields equations of motion for the metric, E^^ = 0, and matter fields, 
Em = 0, via 

6L = Ef6gab + 'Em6iP + d&{cP,6<P). (21) 

For convenience we set E^^ = eE^'' and E^ = eEm- As discussed above we can compute J[^] (defined 
by (0))' define Q[^], for the theory described by (^): it must have the form given in (|lC|): 

J[e] = -2e-Eg-C-€-E^-i;-C + dQ[C] (22) 
(the factor of two between the terms with equations of motion here is purely a matter of convention). 

We can also use the individual Lagrangians and to define the stress-energy tensor T"'', and 
symplectic potentials &g{g,6g) and @m{4',S4')- 

6Lg = E'^'dgab + d@g{g,6g) 

6Lm = Em6iP + €^T''^6gab + d&mi<l),6(l)). (23) 

Clearly E^^ = E'^^ + e^T""^, and up to the ambiguities present in the symplectic potentials, we also 
have = @g + &m- Similarly, if we define the Noether currents for the individual Lagrangians by 

^ @g{g,C^g)-C-Lg, 

Jm[e] = @mi^,C^(p)-C-'Lm, (24) 

then it follows that 

m=^9[i]+Jm[^]- (25) 



Now we impose the structure ( |10D on each of Jg and Jm, in the process defining Qp and Qm, which 
are the Noether charges in the theories arising from these Lagrangians: 

J,[^] = -2e - E'g - ^ + dQg[(] 

3m[^] = -e-Em-^p-^-e-T-C + dQrr,[^]. (26) 
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Finally we substitute ( |2^ ) into the right side of (25) and ( |2^ ) into the left side, obtaining 



-2e-Eg-C-e-E^-i;-C + dQ[C] = -2€-E'g-C-e-T-C-e-E^-ij-C + dQg[C] + dQm[C]. (27) 

All the terms involving equations of motion and stress-energy tensors can be seen to cancel, and the 
resulting identity implies 

Q[^] = QaiC] + Qm[^] + dZ, (28) 

(where Z is some arbitrary covariant one- form). We therefore have a relation (independent of any 
field equations) between the Noether charge, Q, of the full theory given by L, and that of the "pure 
gravity" theory Qg, arising from Jjg. We are now ready to state the identity: 

Lemma 2: Fix L, (the "vacuum" Lagrangian) and L^, (the "matter" Lagrangian) to be diffeo- 
morphism invariant Lagrangians related as given in (^) with the functional dependence shown there. 



Fix a smooth vector field let be defined by (|23|), let Qg[^] be the Noether charge defined by (26) 



for the theory described by L^, and let T"'^ be the stress-energy tensor of the matter fields defined by 
(|23|). Now consider an asymptotically flat, stationary, axisymmetric black hole solution with bifurcate 
killing horizon, in the theory described by L, with stationary killing field (with unit norm at the 
sphere 5"°°, at spatial infinity), and axial killing field 93" (with closed orbits), so that and (p"" Lie 
derive the metric but not necessarily the matter fields. Let the horizon killing field (which vanishes 
on the bifurcation sphere Ti.) be given by x'^ = ~^ ^H^"", where is a constant. Then for 6 a 
perturbation to an arbitrary nearby solution, such that J^** = 0, 

/ 5Qg[^] -^.@g= [ 6Clg[x] -nn [ 5Qg[ip] + I \i-e T'^'dgab - 5(e • T • 0- (29) 



Proof: 

We evaluate the expression ([T8| ) for the theory (p!9|), where the background solution is a black hole 
with the symmetry and structure described above (p9[), demanding that the metric be stationary and 
axisymmetric in the background spacetime, but placing no restrictions on the matter fields. In this 
case the integrand on the left side of (|l^) is generally nonvanishing. Assuming that the field equations 
hold in background for the matter fields, E^ = 0, and that 5^ is ^ solution to the linearised matter 
equations of motion off this background {5 Em = 0), we find the left side of (|^) is 



u>{<i),C(.<j),5(t)) = 5&g{g,C(^g) - C(.®g{g,5g) + 5®rn{(t',C):<j)) - C(.®rn{(t 

= 5{dClm[i] - e • T • e + e • L^) - £50„((/., defy 

= 5{d(^m[i] - e • T • e + e • L^) - e • a!e„(0, 5cP) - d{i ■ &m{<P, S(l>)) 

= d{5Clm[^]-^-&m{<P,Scl)))-6ie-T-^) + ^^-eT-%ab, (30) 



where we used the stationarity of gab in the second line, the expression ( pq ) for in the third, the 
Lie derivative identity C^X = ^ ■ dX + d{S, ■ X) (which holds for an arbitrary form A) in the fourth 
line, and the definition ( |2^ ) of 0m and the stress-energy T"^ in the fifth line. Now also assuming 
Eg = 6Eg = 0, and substituting (|30| ) into the left side of (|l^) yields 

/ di6Qm[C]-C-@n.) + U-eT''''5gab-6ie- T-O = I 5q,[i]-^-&{ct>,5<t>)- ( 5q,[i]-i-@{ct>,5<t>), (31) 
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and so, cancelling the boundary terms 5Qm[C] ~ C ' ©m from both sides (and using (28)) we get 

f U-e T'^'Sgab -6{e - T-0= f SQg[C] - ^ ■ ®g{g, 5g) - f 6Qg[C] - ^ ■ @g{g, 6g). (32) 
2 Js^ Jh 

Now writing in terms of a-nd 99" at the boundary TC (and discarding terms which vanish as a 
result of the vanishing of x"" ^.t 7i, or which vanish because if"" is tangent to 7i) we get 

/ • e T-%ab -5{e - T-0= I 5Qc,[e] - C • ®g{g, Sg) - f 5Q,M + f ^dgM, (33) 
which is what we wished to show. □ 



The identity ( |29D has physical significance when we can interpret the surface integrals appearing there 
as (variations of) the energy, entropy and angular momentum of the black hole. When is this possible? 
If the theory had no matter fields then we could choose to vanish, and the terms involving T""^ in 
( p^ ) would vanish (we could also choose other breakups of L, and we'll return to this shortly). In this 
case we'd have Lg = L, 0^ = 0, and = Q. If in addition there existed a three-form B, (local in 
the dynamical fields, the flat metric r/afe, and its associated derivative, d, at spatial infinity) such that 
at spatial infinity, ^ • 0((/>, Scj)) = ^ • 5B, then (]29| ) can be written 

5£ = ^5S + nH5JH, (34) 

ZTT 

where the varied quantities in (|3j) are defined below, and have well-known physical interpretations 
IQ]. These are (i) The canonical energy of the system, which we define as 

I QK]-C-B, (35) 

(ii) The entropy 5 of the black hole; by taking the functional derivative of the Lagrangian with respect 
to the Riemann tensor (treated as an independent field) we know (setting e^f, to be the binormal to 
the bifurcation sphere) that 

'5/^Q[x] = |^'J5, (36) 

where 

S = -2tt / — eabGcd, (37) 

Jn OUabcd 

and K, is the surface gravity of the background black hole horizon, (iii) The angular momentum of the 
system measured at the black hole, defined by 

Jh = - f Q[^]. (38) 
Jn 



In fact, the angular momentum can be measured either at the black hole horizon or at spatial infinity; 
since the metric is axisymmetric with axial killing field tp"", it can be seen from (0) that J[ip] vanishes, 
when pulled back to a slice to which (p"" is tangent. This ensures (integrating the relation J[ip] = dQ[(p] 
over E) that, for the background solution, 

/ = / QM. (39) 
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In addition, by considering the identity 



cj((/>,50,£^</.) = / 6Q[^]-^-@, (40) 

E JOE 

we see that when Lie derives all the dynamical fields, the left side of this equation vanishes. Since 
is tangent to the two-spheres Ti. and 5"°°, the pullback of the second term on the right side vanishes. 
It follows that 



6 / Q[f]=6 (41) 

Therefore, in spacetimes which have axisymmetric background configurations, the angular momentum 
measured at the black hole is equivalent to the canonical angular momentum J', measured at spatial 
infinity 

J^-f QM, (42) 



s° 

both when ip is an axial killing field, (in the background solution) and for arbitrary solutions which 
are perturbations, 6(j), of the axisymmetric solution. This calculation also shows that the definition 
of is gauge independent, for arbitrary perturbations of an axisymmetric solution. This is because 
6Jh = when we choose 5(j) to be pure gauge, which we see by first setting S(j) = C^(j) for some smooth 
v"", and then replacing 5(f) with a gauge transform 5'(j) which coincides with dcj) in a neighbourhood of 
the bifurcation sphere, but vanishes in a neighbourhood of spatial infinity. Then we have for every Scj) 
(using (|1D), 



SJh = S [ QM = 6' f QM 
= 6' f QM 



= 0. (43) 

So we have that when T"^ vanishes (along with L^), the interpretation of the terms in (|2^ ) is straight- 
forward, and one obtains a formula ( p4[ ) which (bearing in mind the equivalence of J and Jh) is the 
formula (0). 



What if the set of fields ^p is non-empty ? In general, the ambiguity ( |2Q| ) in breaking L into Tjg and L-m, 
stops us from meaningfully interpreting the surface terms in ( [29| ) as perturbations of mass, entropy 
and angular momentum: even if the overall theory is fixed, every choice of generates a different 
relation, with different choices of Q^ etc. We therefore seek more restrictive assumptions under which 
we might successfully identify the surface terms in ( p9| ) . One approach is to fix a particular choice of 
lig and think of it as specifying an independent theory. We assume there exists a form such that 
at spatial infinity, ■'Bg) = ^ ■ @g, and consider the functional Mg defined by 

Mg^ f Qg[^]-C-Bg. (44) 

If we now require that the stress-energy of the matter distribution falls off sufficiently rapidly at spatial 
infinity, such that (near spatial infinity) the metric for any solution of the L-theory approaches a metric 
solution of the L^-theory, and Mg yields the same result on both metrics, then it makes sense to define 
the mass of the system as Mg. We note that if we can also find a form B(i;^) for the full theory, such 
that at spatial infinity ■ B) = 6(j)), then we can also define a canonical energy, for the full 
theory given by (|35|), and in general £ ^ Mg. 
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Therefore, when the stress-energy of the matter distribution falls off sufficiently rapidly, we can inter- 
pret the left side of (p9|) as the variation of the mass of the system. The surface terms on the right side 
of (^) are (variations of) the functionals that would measure the entropy and angular momentum 
of a stationary black hole in the L^-theory. We might therefore be tempted to interpret them as the 
black hole entropy and angular momentum; indeed, since is the Noether charge of the theory, 
we know from Q that one can define a quantity, Sg, by 

Sg = -2Tr / —^e^hecd, (45) 



such that 

5 

IH 



S I Qgix] = ^SSg. (46) 



One might also define a quantity, Jg^, by 

Jn 

Although we made no assumptions about the axisymmetry of the matter fields, we can show, providing 
the support of T"-^ does not intersect some neighbourhood, U, of the bifurcation sphere, that Jgj^ is 
also well-defined (gauge independent) for arbitrary perturbations of the axisymmetric solution. This 
follows by evaluating the left side of (^), using the fact that the calculation ( |30[) also holds when 
is replaced by ip°'. Taking to be tangent to the spatial slice, Eq. (|^) then becomes 

/ 5Clg[^] + ( 6{e.T.^)= f 6(^g[ip] (48) 

Now, as before, let the perturbation in this equation be gauge, 54> = Sep. Then we again can replace 
the perturbation on the right side with an equivalent gauge change, which vanishes outside U, and so 
intersects neither the support of T"'' nor spatial infinity. Then we have the left side of ( |48| ) vanishes, 
and so 

SJgH = - f SQg[cp]=0. (49) 
Therefore Jgj^ is defined for arbitrary perturbations of an axisymmetric solution. 
Now having defined Mg, Sg and Jg^^, we could write out (|29|) in the form 

6Mg 



^5Sg + nnSJg^ + / k • e T'^'Sgab - 5{e ■ T ■ 0, (50) 



where (5 is a perturbation to an arbitrary nearby solution. However, we caution the reader that 
the identification of black hole entropy with Sg in general gives results in conflict with those in 
consider a theory of gravitation with a scalar fleld, for which the matter Lagrangian couples to the 
spacetime curvature, and which displays stationary black hole configurations in which the scalar field 
has sufficiently rapid spatial falloff. We can therefore write out (^0|) and interpret the black hole entropy 
as Sg. From the results of Q we expect the entropy of the black hole to include contributions from the 
scalar field; equation (p0|), however, defines a black hole entropy Sg with only metric contributions, 
with the entropy contribution of the scalar field somehow distributed in the volume integral of its 
stress-energy. These two points of view are contradictory; therefore, while there are clearly special 
cases (for instance, the Einstein-perfect fluid system) in which we can identify Sg as the black hole 
entropy, and terms in the volume integral as (variations of) the matter entropy, in general we regard the 
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notion of the black hole entropy defined by Sg as inappropriate. Clarifying when Sg can be correctly 
interpreted as black hole entropy is the subject of future research. 

We note parenthetically that we can write out an alternative form of (|50|) by replacing the stationary 
killing field in (|l8|) with the horizon killing field x"- (The analysis up to (^) is unchanged except 
for the substitution — > x""-) Then expanding = + ^Hf"" at spatial infinity and on the slice S, 
(but not at 7i) and making the definitions discussed above gives the identity 

6Mg = ^6Sg + QndJg^ ^J^l^'^ ^"''^^'^^ -6{e- T-C)-nH Jj{e- T-ip) (51) 

where Jg^ = — J goo Qg [(/?], is the system angular momentum measured at spatial infinity. Therefore 
the cost we have incurred for the transfer of the angular momentum integral to spatial infinity is the 
appearance of an extra term in the volume integral. 

A relation of the form (|5^) , was first given by Schutz and Sorkin , in the case where was fixed to 
be the Lagrangian for general relativity, was any matter Lagrangian, and there was no black hole 
boundary 7i, for the hypersurface E. The relation stated in is correct, but we comment here on the 
ambiguity of the "Noether operators" used by Schutz and Sorkin to derive it: In its initial definition 
10] the Noether operator for a Lagrangian L and a smooth vector field was defined to be any (not 
necessarily covariant) three form J''^[^] satisfying the relation 

£5L = E£^(/. + d(J'^[e]+e-L), (52) 

for every smooth field vector field This definition leaves 3^[(,] ambiguous by an arbitrary exact 
three-form which is a linear differential operator in Since we know from 1^) that J'^ = dQ[^] 
when the field equations hold, this ambiguity would permit J'^ = as a valid Noether operator 
(which, following Schutz and Sor kin's methods, would yield a correct but trivial relation). On the 
other hand, our definition of the Noether current admits a limited set of ambiguities (stated after 
([Tsl)), which cannot be used to annihilate the Noether charge, and in particular do not change the 
content of the first law. 

Sorkin introduced an augmented definition of the Noether operator in |^], requiring that for a variation 
of the dynamical fields given by 64> = Z^^^, where / is any function, the Noether operator J"^ be 
defined by 

<5L = E/£5(/< + d(/J^'[e] + /^L). (53) 

Providing one can find a J*^' which satisfies this relation, it is easy to see that one cannot add a term 
to J'^ which is both exact and linear in /, for arbitrary /. For a theory with a first order Lagrangian, 
finding such a J'^ is always possible: in a first-order (noncovariant) Lagrangian for Einstein-Maxwell 
theory was used to yield an unambiguous Noether operator. It is not clear, however, that any general 
Lagrangian theory has a first order Lagrangian formulation, so in general, Sorkin's definition may 
not even yield a Noether operator. In contrast, all of our Noether currents J[^] defined above can be 
computed for Lagrangian theories of arbitrary derivative order, and are manifestly covariant, requiring 
no additional background fields (apart from the symmetry field i^") for their definition. For these 
reasons, we feel that whilst our relation ( pO] ) and that in coincide for an Einstein-matter system 
without the black hole, ( |50|) is defined more generally. 

We finally remark that we could have carried out the entire analysis leading up to (p9| ) allowing the 
Lagrangian to depend on a set of stationary axisymmetric fields, Si, including the metric, and the 
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Lag rcLiigicin Lrn to dGpGnd on Si cinci a distinct set of fields, tp^ wliich didn't appear in L^, to obtain 
a relation very similar to (^). The resulting perturbative identity has the terms Qp and &g in (|29|) 
replaced with the Noether charge and symplectic potential in the theory described by (which now 
depends on both the metric and the other matter fields in the set Sj), and the volume term is now 
given by 

f b.eTsM-H^- Ts-s-0 (54) 

where the first term in the volume integral is defined by the variation of L^: 

1 

6Lm = E^^V' + -Tsjsi + d@^{(t), S(t)). (55) 
Giving each Sj field Ui upper, and di lower indices, in the manner 

Si^Si'-^-'-'a^-a,^, (56) 

the second term in the volume integral is defined by 



{e ■ Ts ■ S ■ ^)abc — Ceabc ^[^Sib^...b^. ^ ( Si aj^...a^.6p . . . S,;^ 



(57) 



3 A review of perfect fluids, and three variational formulations. 



In this section we recall the definition, the relevant properties, and three variational principles for a self- 
gravitating perfect fluid: one given by Schutz , (which we use in the Appendix to derive a conserved 
current for perturbations of Einstein-perfect fluid systems), the "axionic vorticity" formulation given 



by Carter [10| for an isentropic perfect fluid (which we use in the next section, to derive a first law), 



and a "convective" approach also described by Carter |10|. Our aim is to gather the results we need 
for the calculations of the following sections; detailed treatments of these variational principles can be 
found in |g. 

From the viewpoint of black hole mechanics, we would like a stationary axisymmetric black hole con- 
figuration to be represented by a Lagrangian theory in which all the fields appearing in the Lagrangian 
(the dynamical fields) are also stationary and axisymmetric. Having stated these formulations, how- 
ever, we will see that they all have fluid configurations in which the physical fields (the fluid four 
velocity, number density, entropy and functions of these fields) are stationary and axisymmetric, but 
in which the dynamical fields possibly share neither of these symmetries. The question as to whether a 
variational principle exists that always represents (physically) stationary axisymmetric configurations 
with dynamical fields that also have these properties is (as far as we are aware) open. 

By a perfect fluid on a fixed spacetime background jl^, [l^ we mean a system described by five scalar 
fields, {n, s, p,p,T), on spacetime and one (unit, timelike) vector field U"", such that p = p{n,s) is a 
fixed function, and the following equations hold on the fields: the first law of thermodynamics, 

dp(n, s) = dn + nTds, (58) 

n 
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and the equations of motion, 

VainU") = 0, and VaT"* = 0, (59) 

where T"^ is defined by 

T"'' = {p + p)U''U^ + pg"'' . (60) 

The fields n, p, T and U°- have physical interpretations as the number density, energy density, en- 
tropy per particle (specific entropy), pressure, temperature, and four velocity of the fluid, respectively. 

We note that (|59|) can be given a useful alternative form, by first defining the specific inertial mass /i: 

II = (61) 

n 

which along with ( [5^ ) implies 

dp = nd/i — riTds. (62) 

By using these relations in the second equation of (|59|) we get (see ) an equivalent pair of equations 
of motion - 

VainU") = 0, and nU^uJab = nTVbS, (63) 
where the fluid vorticity two- form, ujab, is defined as 

0Jab = 2Vla{fiUb]). (64) 

If desired, one can define the entropy per unit volume, 5 (entropy density), by S* = ns. Substituting 
this definition of S into ( [5^ ) and defining the chemical potential, by 

,_p + p-TS 

p = , (65) 

n 

then gives the relation 

dp{n, S) = fi'dn + TdS. (66) 

We now specify three variational formulations for this perfect fiuid, over a fixed spacetime background 
(coupling the theories to gravitation amounts to adding the appropriate metric Lagrangian, which we 
do later). First, we state the "velocity-potential" representation of Schutz here the dynamical 
fields of the fiuid are given by scalars $, a, /?, 0,and a. One now defines a function m which depends 
on these fields via the relation 

= - (Va$ + aVaP + OVaa) (V"$ + aVp + VV) , (67) 

and the fiuid Lagrangian is given by 

Lf = eP{m,a). (68) 
where P{m, a) is some fixed function. One can verify [|, 111 that we recover (m), and also that the 



equations of motion for the fields ^,a, f3,d,a arising from this Lagrangian reduce to (|63|), provided 
one defines the physical fields in terms of the dynamical fields in these equations by: 



p - 


p 


m - 




a - 


-> S 


- 


n 


)m ~ 


-> —nT 


a<7 - 





(69) 
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Conversely, given any configuration of the physical fields {n, p, s,p,T,U°') satisfying ( [Sq ) and (|59|), it 
can be shown (see [^) that there exist functions {P, m) and (non-unique) dynamical fields (<I>, a, /?, 9, a) 
related to the physical fields by (|69|), which satisfy the equations of motion arising from Lagrangian 



Next, Carter's variational formulation [10| for an isentropic perfect fluid, (by which we mean that the 
fluid has an everywhere constant specific entropy s), defines the dynamical fields to be a two- form and 
two scalars, bab and x^- The fiuid Lagrangian is given in terms of these fields by 

L/ = i-riu) - le'''''XbVcX^VdX-]e. (70) 

where the function r(iy) is fixed, and the function ly is defined in term of the potentials bab by the 
relation 

= l{V[abbc]){vH''^^ (71) 



As shown in [10|, if one defines the physical fields as follows: 



r ^ p 

V ^ n 

v{dr/du) — r ^ p, 

3V[Xb] ^ Nabc, (72) 

where the number-density three-form Nabc is given by 

Nabc = neabcdU'^, (73) 

then we recover (|58D, and the field equations for bab and x^ yield the second equation in ( p^ ) in the 
case VaS = 0, as well as the relation 

UJab = 2V[aX^Vb]X'. (74) 



Given relation (|r2D between bab and Nabc, one sees that the first equation in (|63D is satisfied vacuously, 
since it can be rewritten as 

^laNbcd] = 0, (75) 
but the definition of Nabc shows dN = ddh = automatically. 



A third type of variational formulation given by Carter [IC], and treated in more detail by Brown |14], 
(which is the equivalent diffeomorphism invariant version of the formalisms specified by Taub |l^, or 
Hawking and Ellis [^), has dynamical fields for A = 1,2,3. In this formalism one must specify 
two functions- r(u, a), and (t{X), where u is defined in terms of the by 

i^^ ^ 6[iVABc(^)V„X^V6X^VeX^][Afz3EF(X)V"X^V^X^V^X^], (76) 

and Nabc{N) is a fixed three- form on the three-dimensional manifold which has X"^ as coordinate 
fields. The Lagrangian is then given by 

Lf = -er{u,a). (77) 
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The equations resulting from this Lagrangian for the fields X are seen to reduce to the second 



equation in (|6^) after one has set 



r — > /9 

— > n 

cr — > s 

i'{dr / dv)(j — r — > p 

{dr/da)^ ^ nT 

NABc{X)VaX^VbX''VcX^ ^ Nabc, (78) 

where Nabc is defined from (|73|). (This relation between the physical A^afec and the dynamical fields 
also ensures that Nabc is automatically conserved.) The X^ are interpreted as coordinates on a "base 
manifold", obtained by treating the spacetime as a bundle with fibres given by the integral curves 
of the four-velocity. We will not use this formulation for two reasons: firstly, the assignment of the 
entropy, s, as a fixed function of the X^ only allows us to perturb it by diffeomorphisms of the base 
manifold (for this reason we use Schutz's formalism for the calculation in the Appendix). Secondly, it 
is unclear that there are any solutions in which the X^ are globally well-defined axisymmetric fields 
on spacetime (for this reason, in section 4, we use the formulation due to Carter with Lagrangian 



In order to write the first law in form (|2|), only involving surface integrals, we must assume that all 
the dynamical fields are stationary and axisymmetric in the background solution. Now even if a fluid 
configuration has stationary and axisymmetric physical fields (the fluid number density, entropy and 
functions of these fields), the dynamical fields (the fields appearing in the Lagrangian) corresponding 
to these physical fields may not possess these symmetries. Therefore, the requirement of stationarity 
and axisymmetry on the dynamical fields may restrict the choice of background configurations. In 
fact, for Schutz's formulation, we see from the definition of the four velocity ( |69|) that physical fluid 
configurations with an everywhere causal four-velocity (including those which are stationary and ax- 
isymmetric) must include at least one nonstationary dynamical field. There are therefore no physically 
interesting fluid configurations in which all the dynamical fields in this formulation are stationary. 



On the other hand, for Carter's formulation, it is evident that there must be some physically stationary 
fiuid configurations with stationary dynamical fields; (for instance, a static spherically symmetric fluid 
distribution could have the field bab given by b ~ f{r)'^e and = 0, where is the volume element 
on the spheres of symmetry). However, we will see in the next section (in the discussion above (|95[)) 
that a stationary, axisymmetric, circular flow (in a spacetime which also has these symmetries) must 
be vortex-free, if are restricted to be stationary and axisymmetric. That is, the assumption of 
stationarity and axisymmetry on the vorticity potentials restricts the allowed stationary axisym- 
metric configurations a fluid can adopt. We make no attempt here to enumerate the set of physically 
stationary and axisymmetric configurations which also have these symmetries in the dynamical fields 
(or indeed, in the case of black hole spacetimes, to investigate whether this set is non-empty). Rather, 
in the following section we will assume the potentials are stationary and axisymmetric, and write out 
the resulting first law involving only surface terms, looking for any non-trivial modifications arising 
from the fiuid fields. 



We are unaware of a variational formulation for a perfect fiuid which represents all stationary ax- 
isymmetric fiuid configurations with stationary axisymmetric dynamical fields. If it exists, then the 
following argument by Schutz and Sorkin Q shows that certain compactly supported perturbations of 
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the physical fields must correspond to non-compactly supported perturbations of the dynamical fields. 



Since the calculation given in (30) does not depend on the fulfillment of the field equations for gat, it is 
still valid if we consider the fields tp to be the dynamical fields for a perfect fluid over a fixed spacetime 
background, and we let 61^ be a perturbation to a nearby solution of the perfect fluid equations, with 
6gab = 0. Now consider a formulation for a perfect fluid where, for a general configuration in which 
all the physical fields (and the metric of the spacetime background) are stationary, all the dynamical 



fields are also stationary. Then the left side of (30) vanishes, and integrating the right side over a 
spatial slice S, we are left with 

[sie-T-0=[ SQ^[^]-C-@^{<p,5^). (79) 

This implies that for perturbations of the physical fields for which the corresponding perturbations of 
the dynamical fields are compact, we must have 

Jj{e-T-O = 0, (80) 

which, for a perfect fluid, is clearly false for a general stationary background. This implies that if a 
variational formulation is to have dynamical fields which are always stationary when the physical fields 
are stationary, then perturbations of the physical fields which yield a non-zero result on the left side of 
(|7^ ) must correspond to spatially non-compact perturbations of the dynamical fields. This requirement 
rules out the existence of a variational principle in which the physical fields are the dynamical fields 
10]. However, the existence of a variational principle for a perfect fluid in which all configurations with 
stationary and axisymmetric physical fields are represented by dynamical fields with these symmetries 
is still an open question. 



4 First laws of black hole mechanics with perfect fluids 



We now present two forms of the first law of black hole mechanics which incorporate perfect fluids. The 
first form is a special case of the perturbative identity (|5^) , where is the usual Hilbert Lagrangian 
for general relativity, and is any Lagrangian for a perfect fluid. This form of the first law allows 
non-stationary dynamical fields, at the cost of having volume integrals in the interior of the spacetime. 
We then compute a second form of the first law only involving surface integrals for both metric and 
fluid fields, using Carter's variational formulation presented above, and the methods of [Q]. 



4.1 The first law with volume integrals 

We now write out the perturbative relation (|50|), setting hg = 1/16ttR, and to be any perfect 
fluid Lagrangian which allows all possible perturbations of the physical fields of the perfect fluid off 
an arbitrary background. (From the comments below Eq.(^) it is evident that Schutz's variational 
formulation, with Lagrangian ( |6^ ) satisfies this criterion.) As stated in Lemma 2, we assume the metric 
of the background spacetime is asymptotically flat, stationary and axisymmetric with a stationary 
killing field and axial killing field f'^. We also assume the existence of a bifurcate killing horizon, 
with horizon killing field = + ^Hf"", where is the angular velocity of the horizon. 
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In this case (see 0) the term Mg in ( |50| ) can be shown to be the ADM mass, Sg to be 1/4Ah, and 
Jgjj the expression Jh for black hole angular momentum given in (|l]). The terms involving the stress- 
energy tensor have been shown by Bardeen, Carter and Hawking |jl| to reduce to the fluid terms in 
dl]), but for completeness (and to fix the signs for our choice of orientations) we briefly demonstrate 
this fact: in [Q] the four velocity of the fluid with angular velocity Q (which need not be constant) 
was set to be U"" = v'^/{—v ■ v)^/"^ where v'^ = + Q.ip'^. Now using (|60|) , (|6^) and (p6D, (assuming, as 
usual, that we identify the perturbed spacetime such that 5^" = 6(j)'^ = 0), one obtains 

= v'^Siifl'n + TS)Va{-V ■ vY^I^U^Ebpgr + peapqr) " mT\^' e^pgr) 
= {P + pKdiVai-V ■ v)-^/^)U''ebpgr + ^pg^'^SgcdC^eapqr 

-li'{-v ■ v)^/^6{nU''ehpgr) - T{-v ■ vfl^iSU^'e^^,,) 

-{n5fl' + S5T)Ceapqr + V^Sapgr^p - fl5{T"- b^f'' eapqr) 

= eeapqr^T'^%cd + t^' {'V ■ vflH{nU' e^gr,) - Ti-v ■ v)^/H{SU'e,pgr) 

-n5{T\ip''eapqr), (81) 

When all these substitutions are inserted into (^0|), it reduces to 

SM = ^6A + QhSJh - I li'\v\5Nabc + / ^SJabc + I T\v\5Sabc, (82) 

which is identical to (|l|), except that 6 now represents an arbitrary perturbation (not necessarily 
stationary or axisymmetric) of the background. In this sense, (|8^ ) is a generalisation of (||). 



4.2 A (restricted) first law with surface integrals 



In the previous section we observed that the variational formulations we presented were constrained 
in the stationary axisymmetric fluid configurations they could represent, given the requirement that 
their dynamical fields obeyed these symmetries. One might therefore suspect that any form of the first 
law involving only surface integrals could not include non-trivial fluid contributions. Indeed, if we add 



Schutz's Lagrangian (68) to the Lagrangian of an arbitrary metric theory of gravity, and construct a 
first law using the analysis of then we find no additional contributions to this first law from the 
fluid fields, providing the fluid's number density decays sufficiently rapidly at spatial infinity, and does 
not intersect the black hole horizon. It is possible, however, to convert some of the volume integrals in 
dl]) into surface integrals, by choosing Carter's variational formulation (|70|). We do so below, finding 
a first law for an arbitrary metric theory of gravity coupled to an isentropic perfect fluid, in which the 
background configuration for the perfect fiuid as well as the allowed perturbations of the physical fields 
are restricted. (Note that the gravitational contributions to such a first law have been considered in 
detail in Q. We are interested in the fluid contributions.) We finally verify that this first law reduces 
to dl]) when the assumptions made in the two derivations overlap. Our first law is the following result: 

Lemma 3: Let L, given by 



1 



Lg - e{r{u) + ^e'^''^''6„bV[,x+V,]X-), (83) 



2 
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be the Lagrangian for an isentropic perfect fluid coupled to an arbitrary metric theory of gravity, 
where Lg = eLg{gab, Rabcd,^ Rabcd, ■ ■ ■ , iy^)Rahcd), and the perfect fluid formulation, with dynamical 
fields {bah, X^)) is summarised below ([70|). Fix an asymptotically flat black hole solution with bifurcate 
killing horizon, with the spacetime structure and the killing fields described in Lemma (2), with the 
additional assumptions that all the dynamical fields, (not just the metric) in this theory are stationary 
and axisymmetric, and that all the dynamical fields are globally defined. Let dcj) be a perturbation of 
the dynamical fields, from such a solution to an arbitrary nearby solution, with 6^°" = 0. With these 
assumptions the following identity is the first law of black hole mechanics for this system: 



'g = -^SS + QhSJh + IJ-ooSbqr,- IJ-ooSbqr,+ ^gr " / Xqr (84) 

zvr Jn Js°° JH JS°° 



Here we define the mass of the system, Mg, as 



Mg 



f Qg[C]-^-Bg (85) 



and the entropy, S, and angular momentum, J'h, of the system by 

S = -211 [ f^'^ eabGcd 
Jn OKabcd 



Jh ^ - Q,[vp], (86) 
Jn 

where k is the surface gravity of the black hole, the two- form Qg[^] was defined in ([2^), and the 
three- form Bg is such that, at spatial infinity, ■ B^) = ^ • 0^, with 0^ given by ([2^). Finally, the 
two-form Xqr is defined by 

Xqr ^ 2ebp[q[6{flUr^) - VrlX^X^ + V^lX+^X"], (87) 

Proof: 

The first law of black hole mechanics in is essentially given by the right side of (|l8|) , when the left 
side vanishes because of the assumed symmetries of the background fields. We therefore compute the 
quantities appearing in the right side of (18): Varying the dynamical fields in L (and performing the 
substitutions (72) where applicable) yields the equations of motion and the symplectic potential 0: 

f6Lg 1 ab\ x„ ,^^AT . ^abcdi 



\ogab 2 y 6 

+e[V,{^N-'') - ]e-''^L0,a]6bab + d&, 
In 4 



with the stress-energy tensor 
and the symplectic potential 



rj^ab ^ J^j^a^^^bcd _ pgab^ 

2n 



@pqr{^, 6cl>) = @gpg,{g, 6g) - ^N^^'dbbceapqr + ^eapqrbcde'"''^ {V bX^^ ^X' " ^bX'Sx^)- (90) 



It can be verified that the equations of motion for the fluid fields reduce to (|63|) using the definitions 
(|7|,|73D. The stress -energy tensor ( |89[ ) is also seen to reduce to the usual form (^) by expanding its 
first term: 

^j^a^^j^bcd ^ ^^jje^a^^^^bcdf^jj^ ^ ^^^^ab ^ ^ajjby 
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The Noether current associated to is 

Therefore, the integrand on the right side of ( p^ evaluates to 

Ve(VpX~<5x+ - Vpx+<5x-) + (92) 



where we define the two- form Xqr by (|87|), and we used the identification Nabc = ^abcd^U to obtain 
the second hne of (lol). 



When the background solution is a black hole with the structure and symmetries specified in the 
statement of the Lemma, the fourth term in the second equation of ( |92| ) vanishes because the dynamical 
fields are stationary: • = 0. Now given the definition of vorticity (|6^ ) and its relation to the 
potentials (|7^), it is evident that (locally) there exists some function / such that nUa can be rewritten 

/xC/a = Va/ + x+V,x"- (93) 

Let t be a function such that £,'^dta = 1. Then the requirements that the four- velocity be causal, 
stationary and axisymmetric, along with the assumed stationarity and axisymmetry of force / to 
be a sum of terms, one of which is strictly linear in t (we define the constant of proportionality to be 
— /Ucxd). For the same reason the (^-dependence of / must be also linear, but this dependence can be 
ruled out because the occurrence of such a term would force to be acausal near spatial infinity. We 
therefore have that the form of / is 

f = -f^oot + g. (94) 

where ^ • = ■ Vg = 0. Therefore we see that the assumption of stationarity and axisymmetry 
on the dynamical fields (taking the four-velocity to be everywhere causal) has restricted us to a very 
narrow range of allowed background four-velocities; for instance, we must have '^"'Ua = 0. Moreover, 
when the vacuum theory is general relativity, with the flow assumed to be circular (tangent to the 
^ — (/9 subspaces), there is only one possible solution: for this theory the subspaces orthogonal to 
and 93" are integrable, and the resulting submanifolds can be endowed with coordinates (x^, x^), such 
that the metric is "block diagonal" with no "cross-terms" between the subspace spanned by ip°' and 
its orthogonal complement (see Chapter 7 of [^]). Now the assumption of circular flow forces g = Q 
and x^dx~ = 0, leaving us with only 

liUa = -Hoodta. (95) 
In any case, using just the form of / in (p^), we see 

eVC/a = -eVoodta = -/Uoo, (96) 

and the boundary term (|9^ ) reduces to 

'^Q,r m - e@pgr = SQg^^ ^ - e&9pgr + /"-^V + ^gr, (97) 
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We now assume the existence of a form Bp such that at spatial infinity ^ • 0^ = • Bg), and write out 



the first law of black hole mechanics by substituting (97) into the surface integrals on the right side 
of dl^), observing that the left side of (|l8|) vanishes due to the symmetries assumed on the dynamical 
fields. If we expand = ~ ^Hf"" at the bifurcation sphere for the first two terms of (|97|), then we 
obtain (|8^) which is what we wished to show. □ 



The results of |^ predicted that the first law (|84| ) would only contain surface integrals, and we see this is 
indeed the case. Note, however, that the assumptions made about the symmetry of the dynamical fields 
restricted the allowed background fluid configurations for the fluid fields. Moreover, by perturbing the 
local form of /xC/q in ( p3| ) we see that the restriction to stationary and axisymmetric i^i background 
also prevents us from achieving all possible perturbations of ^Ua-, by perturbing only the dynamical 
fields bah and x^- Finally both the background and the perturbed configurations must be restricted 
such that the integral /^oo Xq^ converges. (This, along with the following result relating this term to 
the fluid angular momentum will guarantee the convergence of the corresponding boundary term at 
the bifurcation sphere). 

We finally show that (|^) reduces to (^ when the assumptions made in the two derivations overlap. 
From our discussion in the last section we know that Mg , S and Jh reduce to their values for general 
relativity given in (^), when Lg = (l/167r)i?. We start by considering the fluid contribution in our 
first law (|8^ ) from the integral 

Joo JH JT. 

M^Npqr, (98) 



where the last line follows because the fluid flow in |Q] is assumed to be tangent to the subspaces 
spanned by and 99": so taking the velocity to be U"" = v"" /\v\ where v"" = + Vtip'^, we see from 
the discussion above (p^) that ji = —fiU ■ U = fiooV ■ dt/\v\ = fioo/\v\, and so fi\v\ = /Uqo- Our first law 
now takes the form 

SM = ^6A + nH6JH - [ MSNabc+ [ Xgr- [ Xgr (99) 



We now concentrate on the original form of the first law in ([T|) and show that it agrees with ( |99|) . By 
repeating the calculation (^Tj) using the relation ( [5^) instead of (|66|) along with the assumption 6s = 
(as befits an isentropic fluid), we find the form of (||) for an isentropic fluid - 

SM = ^6A + Qh6Jh - f MdNabc + / ^SJabc- (100) 

Next, we demonstrate that the pullback to S of the angular momentum density given in ( |100 ) reduces 
to the exterior derivative of the two form Xg^ defined in ( p^ ) , given the assumption that the dynamical 
fields are stationary and axisymmetric, i.e., 

^SJpgr = -{dX)pgr, (101) 

where both sides are assumed pulled back to S. To do this we compute the exterior derivative of (|87|), 
finding 

{dX)pgr = 3eiV,[p,(5(/it/,]) - V,]X-<^X+ + V,]X+<^X"), (102) 
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where we have assumed C^hab = 0. Pulhng this form back to S by contracting with l/Ge^^'^^ns (where 
ng is the unit normal to S) yields 

dX = -^e{2nn,6{fiUr)uH'^ + 2e[^C/"lne(V,x"'^X+ - V^x+^x")), (103) 

where is the volume form induced on S: ^e^cd = n°'^abcd- Now using the axisymmetry of the 
X^, and writing U"" as = v"" /\v\ with angular velocity 17 as given above in (99), we have (using 
5^'' = Scf)" = 0) 

dX = -^e2nn^Q5{iJLUr)ip^^C^ /\v\ 

= -n6J (104) 



where J is the pullback of Jabc to S. Therefore (^) now matches ( p. 00] ) and so the first law (|8^) now 
agrees with the first law given in (||). 



Acknowledgements 



This research was supported in part by NSF grant PHY-95- 14726 to the University of Chicago. This 
work is presented as a thesis to the Department of Physics, The University of Chicago, in partial 
fulfillment of the requirements for the Ph.D. degree. I am very grateful for Bob Wald's expert guidance 
throughout this project. I also wish to thank Bob Geroch, Eanna Flanagan and Shyan-Ming Perng 
for helpful discussions. 



Appendix: The Chandrasekhar-Ferrari conserved current 



The symplectic form uj{((>, 6i4>, ^20) defined in (|16|) is closed when 5i^2(t> satisfy the linearised equations. 
Its dual u)'^{(t), Sicj), S24>), defined by uJabc = ^'^^dabc, is therefore a covariantly conserved current for the 
Einstein-perfect fluid system. Chandrasekhar and Ferrari |]l3[ have, from flrst principles, also derived 
a conserved current, £°'{(j), S4>), for the Einstein-perfect fluid system. Their current is quadratic in the 
(complex) perturbations 5(f), and is restricted to the case where is a static axisymmetric solution, and 
6(j) is a "polar" (even parity) perturbation with harmonic time dependence (we will define this below). 
We now show the equivalence of the uj°'{4>, S(j), dtj)*) and for the Einstein-perfect fiuid system. This 
calculation is the analogue for the Einstein-perfect fiuid system of the calculation by Burnett and 



Wald 1 12] for the Einstein-Maxwell system. 



We start by choosing the Lagrangian for the Einstein-perfect fiuid system to be 

l-pqrs = epqrs[-^R + P{'m,a)], (105) 



where we have set the constant in front of the Ricci scalar to give the field equations in [13|, and used 



Schutz's velocity-potential representation, with Lagrangian (pq). The symplectic potential, 0, arising 
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from this Lagrangian is (after substituting (|69| ) where appUcable), 

©pgr = -\eapgr{yVb " V^t) " n^J'^e.pg, (5«> + a6/3 + 96s), (106) 

where -jab = (^ffafc and 7 = g'^^jab- The resulting presymplectic form is (from (p^)) 

a;p,,(0, <52</.) = -^eapg.[(72^'' - /S) V^Ticd " (272'='' - 72/'')Ve7i/ + 72'^''V,7i] 
-<J2(eapgrnC/")(5i«> + a6i/3 + eJis) 

-eapgrnC/"(<52a5i/3 + 6296is) -(1^2). (107) 



This form is dual to a generally conserved current: it can be shown ||ll|] that for a;" defined above, we 
have (for perturbations and satisfying the linearised field equations), 

VaUj" = 0. (108) 



We now relate this conserved current to the current presented in |lj], by fixing a coordinate system 
with derivative operator da, and writing the volume element e in terms of the coordinate volume 
element e of this system - 

f^abcd = V-geabcd (109) 

- then the vector field w"" defined by Upgr = w^'Gapqr is conserved in the sense daw'^ = 0. If we follow 
Chandrasekhar and Ferrari [^] and specialise to the case where the background spacetime is static 
(with static killing field t") and axisymmetric (with axial killing field (/?"), and the perturbations are 
time and angle-dependent only "harmonically" (that is, there are constants a and lo such that 

CipSrj = iu)5rj, (110) 



for all the dynamical fields rf) then (following ||T2[) it's easy to see that for complex 5(f), w {(j),5(f),5<j)*) 
and w'^{(j),6(j),S(l)*) are independently conserved: dtw^ + d^w'^ = 0. We can therefore restrict our 
attention to the vector components {w'^,w^). Moreover, ( |11[1| ) allows us to substitute the variations of 
the fluid potentials 6{^,f3,s) for variations of their time derivatives: we do this and (recalling (|69|)) 
find 



w'' = w;" - S2i^/^nU'')—6l (fit^Ub) - ^/^—nU''[62 {t ■ Va)5i {t ■ V/3) + <52 (t • Ve)5i (t ■ Vs)] -(1^2), 

fill) 



— / 10 I * / — / I t 1 1 I I / i-i I 11^^^ 

"gr 

where we labelled the contribution from the first two lines of ( |107] ) by Wg^. 

Our aim is now to show the equality of {w'^[(t),5(t),5(l)*),w'^[(j),5(t),5(l)*)) and {£'^,£'^). To do this we 
first specialise the background and perturbations in w"" to those used by Chandrasekhar and Ferrari. 
In the coordinates given in |^] the metric is written 

gab = -e^^dtadtb + e^^dy^adifb + e^^^dxldxl + e^^^dxldxl (112) 

and the nonvanishing (polar) metric perturbations are taken to be 

= -2e2-Jz. 
71^^ = 262^^,5^ 



7122 = 2e^''^5fi2 

7133 = 2e2^MM3, (113) 
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We then set j2ab = Tiafe' where the perturbed functions, 5z^, etc. are complex, but the unperturbed 
functions are real. 

A direct substitution of these perturbations into Wg^ yields the result already known from fl^ ; 

+/^3,2<5(/U + ll^ySins - /^2)] - C.C. (114) 



We now turn to the fluid contributions wj to the conserved current, defined by w'^ = w"" — Wg. We set 
the four- velocity of the background to be = e~^t" and (following |jl3[) and denote the perturbations 
of the orthonormal frame components of U"' by ia^": 5Ua = icr^a- We then find the '2'-component of 
tUj given by 



-dxl — [52{e''+^+^'^+^'■'nU'')5l{neUc)] -(1^2) 
ia 

= ^e'^''+^+^^n62{U^){Sifi + fidiu - /u^iC/q) -(1^2) 

= -e2^+'^+'^3(n5/i + n^(5z^)e; -C.C, (115) 

where we have set 61 = 6 and 62 = 6*, and used the result (see |]l3|) that 617^ = 0. We can also put 
n5fj, = 6p + nT6s, and bearing in mind that 5s must also have harmonic time dependence, we can 
write 

ianT6s = nTt • V6s 

= nTe^U ■ V5s 

= nTe^iSiU ■ Vs) - 5{U) ■ Vs). (116) 



Referring to ( |63| ) we see that the first term on the right side of ( [L16D vanishes whenever the perturbation 
satisfies the linearised equations. Since the background is vortex-free, we see that the second term also 
vanishes as a consequence of (|6^). Adding the resulting fluid contribution to the gravitational terms 



(114) yields 



= _ie'^+V'-M2+M3[5^25(V' + /"3)* 
50,2^(1^ + ;U3)* + 5^l^,25{v + i^Y 

_g2i.+^+M3 + + p)du)Q - C.C. (117) 
Now using the appropriate linear combinations of the linearised Einstein constraint - 

5(V'+/i3),2 + 0,2<5(V' - /i2) + /^3,2<5(^3 " M2) " i^,2<J(V' + ^^3) = 26^+^^ {p + p)^ (118) 
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- to replace the third, fourth and fifth lines of (117) we get 



+2e'^+'^2 ((p + p)5(V + /i3 - ^2)* - 6p*)i2} - c.c. (119) 
where, we define [j4, = A^iA* — AA\. This is seen to agree (up to an overall constant) with S'^ 



of the conserved current in [13|. A similar calculation for w yields £ (which is obtained from £ 
by interchanging 2 <-> 3), and so we find {tv^^w^) = {£'^,£^), and our symplectic current w"' for the 
Einstein-perfect fluid system agrees with the Chandrasekhar-Ferrari current for this system. 



We make two final comments. Firstly, from the comment following Eq.(6£), we know that every con- 
figuration of the physical fields of a perfect fluid has a corresponding equivalence class of configurations 
of the dynamical fields, and as a consequence, every perturbation of the physical fields has a corre- 
sponding perturbation of the dynamical fields. Now, two distinct perturbations of the physical fields 
off the same background (physical field) configuration will each select a corresponding perturbation 
of the dynamical fields. The background dynamical field configuration for each of these perturbations 
will certainly lie within the equivalence class corresponding to the given background physical field 
configuration: however, in general, these background dynamical field configurations will be distinct 
elements of this equivalence class. In using symplectic methods to derive £"" we have implicitly re- 
stricted ourselves to those pairs of perturbations of the physical fields where the corresponding pairs 
of dynamical field perturbations {Si(p,624>) have identical background configurations. In fact, as we 
have seen above, the resulting conserved current agrees with the Chandrasekhar - Ferrari current for 
all pairs of perturbations of the physical fields, not just those restricted in this way. 



Secondly, we notice from (111) that as long as the U"" of the background solution lies in a plane 
tangent to the subspace spanned by and 93", the last term in ( [111[ ) vanishes for the components of 
interest. This in turn yields a conserved current {w'^,w^) which only depends on perturbations of the 
physical fields, without the explicit appearance of the fiuid potentials, for any stationary background 
configuration in which the fiuid velocity is tangent to the t — ip subspaces. Of course, we know that 
iv"" is a conserved current off any background; this observation suggests only that a current similar 
in style to that presented by Chandrasekhar and Ferrari also exists for a background with a fluid in 
circular motion, as well as the static case considered in |13|. 
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